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INTERACTION OF THE MAGNET0R0TATI0NAL INSTABILITY WITH 
HYDRODYNAMIC TURBULENCE IN ACCRETION DISKS 

Jared C. Workman 1-2 and Philip J. Armitage 1,2 
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ABSTRACT 

Accretion disks in which angular momentum transport is dominated by the magnetorotational in- 
stability (MRI) can also possess additional, purely hydrodynamic, drivers of turbulence. Even when 
the hydrodynamic processes, on their own, generate negligible levels of transport, they may still affect 
the evolution of the disk via their influence on the MRI. Here, we study the interaction between the 
MRI and hydrodynamic turbulence using local MRI simulations that include hydrodynamic forcing. 
As expected, we find that hydrodynamic forcing is generally negligible if it yields a saturated kinetic 
energy density that is small compared to the value generated by the MRI. For stronger hydrodynamic 
forcing levels, we find that hydrodynamic turbulence modifies transport, with the effect varying de- 
pending upon the spatial scale of hydrodynamic driving. Large scale forcing boosts transport by an 
amount that is approximately linear in the forcing strength, and leaves the character of the MRI (for 
example the ratio between Maxwell and Reynolds stresses) unchanged, up to the point at which the 
forced turbulence is an order of magnitude stronger than that generated by the MRI. Low amplitude 
small scale forcing may modestly suppress the MRI. We conclude that the impact of hydrodynamic 
turbulence on the MRI is generically ignorable in cases, such as convection, where the additional 
turbulence arises due to the accretion energy liberated by the MRI itself. Hydrodynamic turbulence 
may affect (and either enhance or suppress) the MRI if it is both strong, and driven by independent 
mechanisms such as self-gravity, supernovae, or solid-gas interactions in multiphase protoplanetary 
disks. 

Subject headings: accretion, accretion disks — hydrodynamics — instabilities — MHD — turbulence 



> 

o 

\o 

o 

oo 
o 



X 



1. INTRODUCTION 

The magnetorotatio nal (M R I) or Balbus -H awley 
instability (IVelikovl 119591: IChadrasekhail 119611 : 
iBalbus fc Hawlevl ll99lT 119981 ) underpins the most 
important - - and possibly only - source of out- 
ward angular momentum transport in a wide class 
of well-ionized accretion disks. The MRI destabilizes 
disk flows in which di! 2 /dr < 0, and leads to a state 
of sustained magnetohydrodynamic (MHD) turbu- 
lence that transports angular m omentum outward 
(lHawlev. Gammie fc Balbusl [19 95; Bra ndenburg et al. 
19951:1 Arm itage 1998: Hawlev 20001 iPapaloizou fc Nelson 
2003t lHirose. Krolik fc Stondl2006t) . The majority of the 
transport is mediated by Maxwell rather than Reynolds 
stresses. Following convention, the efficiency of angular 
momentum transport within disks is measured via an 
equivalent Shakura-Sunyaev (1973) a parameter, which 
can be expressed in terms of the fluctuating velocity and 
magnetic fields as, 



aP = pv r Sv,p + 



Mo 



(1) 



where Pq is the thermal pressure. The first term 
on the right-hand side of this equation represents the 
Reynolds (or fluid) stress and the second term repre- 
sents the Maxwell (or magnetic) stress. There is no 
strict equivalence between the evolution of MHD turbu- 
lent disks and models that assume an a shear viscosity 

1 JILA, Campus Box 440, University of Colorado, Boulder CO 
80309; workmanj@colorado.edu, pja@jilaul.colorado.edu 

2 Department of Astrophysical and Planetary Sciences, Univer- 
sity of Colorado, Boulder CO 80309 



(|Balbus fc Hawlevl fl998t iPessah. Chan fc Psaltidl2008T) . 
but for our purposes a is a convenient measure of the 
efficiency of angular momentum transport. 

The fact that the MRI dominates the transport of 
angular momentum within accretion flows does not, of 
course, imply that other sources of turbulence do not 
exist within disks. The most striking example occurs 
in galactic disks, in which turbulence can be driven 
by thermal instabilities and superno va explosions in re- 
gions that are unstable to the MRI (Pio ntek fc Ostrikerl 
2005). However, even in "normal" accretion disks around 
stars or compact objects there are numerous possibili- 
ties, including self-gravity i n sufficiently ma ssive disks 
(jToomrdll96l . convection (jStone fc Balbusl 1 1996f ). and 
Kelvin-Helmholtz instabilities excited by the interaction 
between g aseous and solid components of protop lane- 
tary disks ijCuzzi. Dobrovolskis fc Champnevlll993h . Ad- 
ditional fluid motions can also be driven within disks 
due to the gravit ational influence of embe dded planets 
dBate et"aT1 120021 ) or binary companions (jSpruit et alJ 
[19870, though the wave-like fluid motions induced by 
these sources are different from those initiated by turbu- 
lence. Some of these processes may drive fluid motions 
within the disk whose kinetic energy density is compara- 
ble to that resulting from the MRI. In particular, a self- 
gravitating disk can remain stable agains t fragmentation 
while generating an effective a ~ 10 _1 (Ga mmiel 120011 : 
iRice. Lodato fc Armitagell2005[ ) — as large or larger than 
that produced by the MRI. Since the fluid motions result- 
ing from hydrodynamic processes are uncorrelated with 
the MRI, it is reasonable to suspect that there could be 
non-trivial interactions between the MRI and other tur- 
bulent processes in disks within which both are operating 
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simultaneously. Indeed, iFromag et al.l (|2004f) found that 
in a self-gravitating, MRI-turbulent disk, the strength of 
the angular momentum transport from the self-gravity 
was both weaker, and had a different time-dependence, 
when compared to a disk in which self-gravity alone was 
at work. This re s ult is striking, since the disk simulated 
bv lFromag et al.l (2004) was dominated by low-order (az- 
imuthal wavenumber m = 2) self-gravitating structure 
whose scale was much larger than the most unstable MRI 
scales. 

In detail, different drivers of hydrodynamic turbulence 
may influence the MRI in a unique manner, requiring a 
case-by-case study (self-gravity, for example, is a special 
case since it does yield outward transport of angular mo- 
mentum, whereas most other hydrodynamic mechanisms 
yield negligible or even inward transport). We do not ad- 
dress such subtleties here, but rather study how generic 
driven hydrodynamic turbulence of specified strength in- 
teracts with the MRI. One's expectation is obviously that 
hydrodynamic turbulence that is weak (say, in terms of 
the saturated value of the kinetic energy density) com- 
pared to that driven by the MRI ought to leave the MRI 
unscathed, with significant interaction developing when 
the two sources of turbulence are of comparable strength. 
Very strong turbulence can amplify magnetic fields in the 
disk independent of the MRI, though this may not neces- 
sarily be accompanied by angular momentum transport. 
Our goal in this paper is to test such order of magnitude 
intuition. 

The focus of this paper is on small-scales, which can be 
captured most effectively using local shearing-box sim- 
ulations. Such simulations have significant limitations 
that must be borne in mind. In particular, the strength 
of angular momentum transport in zero-net flux sim- 
ulations with purely numerical viscosity and magnetic 
diffusivi ty has a marked dependence on numerical res- 
olutio n ([Gardiner fe Stonel l2005: Fro mang fc Papaloizoul 
I2007T) - Moreover, when physical values for the trans- 
port coefficients are used (as in this paper) the strength 
of turbulence depends upon the magnetic Prandtl num- 
ber Pm = v/ri as well as on the Reynolds number 
(|Lesur fc Longarettil [20071 iFromang et al.ll2007h . What 
this means for real disks — which when highly ion- 
ized have values of the diffusivity r\ and viscosity v 
that are much smaller than can currently be simulated 
- is unclear, but an obvious implication for numer- 
ical experiments is that the absolute value of a de- 
rived from shearing-box simulations must be viewed 
with caution. For the time being, constraints derived 
from modeling of observed sy stems may be more reliable 
([King. Pringle fc Livid [20CV71) . For this reason, our focus 
here is on how the strength of angular momentum trans- 
port varies in the presence of additional hydrodynamic 
turbulence, rather than its absolute value. 

The plan of this paper is as follows. In §2 we describe 
the set-up for our runs, which make use of the PENCIL 
MHD code previou sly employed for both disk and turbu- 
lence calculations ([e.g. Haugen. Brandenburg fc Doblerl 
2004). In §3 we present results, which concentrate on the 
influence of different levels of hydrodynamic turbulence 
on the saturation level, structure, and angular momen- 
tum transport efficiency of the MRI. These results are 
summarized and discussed in §4. 



2. METHODS 

We use the Pencil Code 1 to solve the dynamical 
equations used to describe a local, magnetized patch of 
an accretion disk in the co-rotating frame. The Pen- 
cil Code solves the equations of compressible MHD in 
non-conservation form using sixth order spatial deriva- 
tives and third order temporal derivatives. The Pen- 
cil Code has been used previously to study the MRI 
(|Brandenburg et aTlll995l ), and has been shown to yield 
similar results to calculatio ns performed using ZEUS 
and similar algorithms (Hawlev. Gammie fc Balbusl 



ll995 HBalbus fc Hawlevl ll998fh One important difference 
between Pencil and ZEUS is that Pencil requires the 
use of explicit diffusive terms to maintain stability. Pen- 
cil is a fixed grid, Eulerian code, that is parallelized us- 
ing MPI. 

2.1. Equations 

We follow the motion of a magnetized parcel of isother- 
mal gas in the co-rotating frame. A local, Cartesian 
co-ordinate system is used, in which x represents the 
radial direction, y the azimuthal direction, and z the 
vertical direction. The equations of motion describing 
the system are derived by solving for deviations in the 
flow velocity u from the Keplerian shear flow U y (x) — 
— 3/2)£lx, where fHs the local background rotation r ate 
Brandenburg et al.l[l995HWisdom fc Tremaindll988l ). 

Using the framework described above, the momentum 
equation is given by 



du 

~dt 



-4 0) 

VP 

P 



du 
dy 
J x B 



P 



(u ■ Vu) + 2Qu y 
F + is h V 6 u 



1 



-flu, 



(2) 



where u is the velocity vector, P is the pressure, p is 
the density, J is the current, B is the magnetic field, and 
Uy is the background Keplerian shear. The terms in 
equation ([2|) correspond to: advection due to the back- 
ground flow, advection of the perturbed velocity field, 
Coriolis and shear effects, pressure gradients, the Lorentz 
force, our forcing term (to be described in the next sec- 
tion) and viscous terms. We neglect the effects of verti- 
cal gravity both for the sake of computational economy, 
and because it is likely irrelevant for this study focusing 
on the interplay between the MRI and forced hydrody- 
namic turbulence. The viscous term ^/jV 6 u describes 
so-called hyperviscous dissipation, specified by a param- 
eter Vh- We use hyperviscosity to limit the effects of 
dissipation to the smallest scales and thereby extend our 
inertial range. For reviews of the effects of hypervis- 
cosity on MHD sim ulations and a derivation of th e hy- 
perviscous term see iBrandenburg fc Sarsonl (|2002h and 
Uohansen fc Kkhrl (120051 ). 

The continuity equation is, 



dp 

dt 



AO) 



dp 
dy 



V-(pu) 



(3) 



where the first term comes from advection due to the 
background flow and the second term refers to the stan- 

1 The Pencil Code and its documentation are available at 
http://www.nordita.org/software/pencil-codc/ 



The MRI in stirred disks 



3 



dard mass flux. We use upwinding rather than diffusive 
terms to stabilize the mass conservation equation. 

The final equation describes the evolution of the mag- 
netic vector potential. Pencil solves for the vector po- 
tential A rather than the magnetic field B = V x A, 
which ensures a solenoidal (divergence free) magnetic 
field. The equation for the evolution of the vector po- 
tential is given by 

¥ = ^ ) | + ^ + ux(VxA) + %v6A ' (4) 

where the first term describes advection of the poten- 
tial due to the background flow, the second term de- 
scribes magnetic stretching due to shear, the third term 
describes the standard electromotive force, and the last 
term is the hyperdiffusive resistivity. 

The final equation needed to close the system is the 
energy equation. We use an isothermal equation of state 
so this is simply P — pc 2 s , where c s is the sound speed. 

2.2. Forcing Term 

Our goal in applying hydrodynamic forcing is to mimic 
the effects of hydrodynamic turbulence that is gener- 
ated within the disk by mechanisms that are indepen- 
dent of the MRI. There are many possible approaches to 
accomplishing this, that vary depending upon whether 
the forcing is applied throughout the volume or only at 
the boundaries, and in the temporal and spatial scales 
involved. The "best" approach evidently depends upon 
the physics one seeks to mimic - forcing that is gravita- 
tionally driven, for example, would logically be written as 
the gradient of a potential function. For most of our runs 
we use a forcing term that injects energy throughout the 
volume at large scales (physically, this would be roughly 
the disk scale height). In the hydrodynamic regime this 
energy then cascades down to smaller scales where it is 
eventually dissipated (by viscosity at very small physical 
scales in a real disk, by hyperviscous dissipation at the 
grid scale in the simulations. We have also considered 
the effect of forcing the fluid on small spatial scales. In 
either case we adopt the forcing function introduced in 
Haugen et. al. (2004) which is given by 

F(x,t) = f k jV e [ ik M- x+i< ^. (5) 

In this function the wavevector k(t) and <p(t)e[— tt, it] are 
chosen randomly at every time step guaranteeing that 
the forcing due to the function is uncorrelated in time. In 
order to ensure that the energy injection rate is indepen- 
dent of the time step is it normalized (using dimensional 
arguments) by setting, 

N=f Q C 3 J 2 y/\k\/Si, (6) 

where the scalar /o defines the amplitude of the resulting 
forcing. We discuss how we set /o in the next Section. At 
each time step a random wavevector k is selected from 
a certain range in k space along with a unit vector e. 
These vectors are then used in the construction of the 
argument of the exponent as well as in the creation of 



which describes a nonhelical, transverse wave. 



2.3. Setup 

We choose a box sized (2tt) 3 with a resolution of 
160 3 . The boundary conditions in the y and z direc- 
tions are periodic and those in the x direction are shear- 
in g. Shearing sheet boundary condi tions are described 
in lHawlev. Gammie &: Balbusl {1995). We initialize the 
system with a density of p — 1, a sound speed c s = 1, 
and an angular rotation rate of fl = 0.2. This scaling 
returns a box size of 2ttH/5. We include hyperviscosity 
to stabilize the scheme at the grid scale. The coefficients 
for the fluid viscosity and the magnetic resistivity are set 
equal to each other yielding a (hyperviscous / hyperresis- 
tive) Prandtl number (vh/Vh) of 1 for most simulations. 
Our hyperviscous and hyperresistive Reynolds numbers 
for grid scale shocks are 75. One technical point to be 
borne in mind when comparing our results with other 
MRI simulations is that we use a cubical domain rather 
than one elongated in the azimuthal direction. A cubi- 
cal domain maximizes the range between the dissipation 
scale and the forcing scale, which necessarily has power 
on the smallest dimension of the domain. The uniform 
resolution in the three spatial directions ensures that 
hyperviscous and hyperresistive dissipation act isotropi- 
cally. We do not see any significant differences between 
our results and prior simulations that we could ascribe 
to the differing simulation domains. 

For our first simulation we initiate the MRI by set- 
ting up a magnetic field in the z direction that has the 
form Ao cos(fcir)z where Aq = 0.05. This yields a plasma 
(3 (defined as the ratio of the gas pressure to the mag- 
netic pressure) of [5 = 800 at the peak of the mag- 
netic field. The instability is seeded by perturbing the 
velocity field with random gaussian fluctuations whose 
magnitude is of the order of 10~ 3 c s . With this setup 
the most unstable wavelength for the MRI, given by 
A c = 2irV a jV3n (jBalbus fc Hawlevl Il99lh . is resolved 
across 23 grid points. We do not use any forcing for the 
initial run which is intended to allow the MRI to develop 
into a saturated state. After evolving the simulation to 
20 orbits we stop it at this point and create a restart 
dump. From here we restart the MRI only simulation 
and let it run for ^130 more orbits. We call this sim- 
ulation 0/o and it is our 'fiducial' run with which we 
compare all subsequent runs. 

2.4. Large Scale Forcing 

Once the MRI simulation is completed we begin sev- 
eral sets of hydrodynamic simulations with equivalent 
parameters but no magnetic fields. To these simulations 
we add the forcing function given by equation ([5]). We 
force the simulations around k = 1.5 to simulate the ef- 
fects of generic hydrodynamic turbulence being input at 
scales slightly smaller than our simulation domain. We 
force slightly interior to the boundaries to avoid issues 
with the boundary conditions. The forcing injects en- 
ergy into the box at large scales, which then cascades 
down in the normal hydrodynamic manner. We run hy- 
drodynamic simulations with varying values for /o until 
the time averaged level of kinetic energy in the hydrody- 
namic simulations, due to the forcing function, approxi- 
mately matches the time averaged level of kinetic energy 
in our MRI only run 0/o- We label our hydrodynamic 
only run as HL and allow it to run for ^150 orbits. 
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Having determining the appropriate value for /o we 
initialize six runs starting from the fully saturated MM 
state we saved at 20 orbits in the MM only case. We im- 
plement forcing from this point onward in each of the 
simulations with values of [~/ Q , l/ , 2/ ,4/ , 8/ , 16/ ] 
and run the simulations for approximately 130 orbits 
each. We label each run by the value of /o; HL is the 
purely hydrodynamical (large scale) run, 0/o is the fidu- 
cial MM only run, .5/o corresponds to the run with half 
strength forcing, and so on. The details of each run are 
summarized in Table [1] 

In addition to the simulations described above we also 
ran four others with some of the underlying physics mod- 
ified. The goal of these runs was to test how robust our 
results are to plausible changes in the forcing function 
(which is essentially arbitrary) and plasma microphysics. 
In two simulations we changed the time scale over which 
the forcing function switches to new wave vectors. The 
function is still defined and normalized as described by 
Equation ([5]) except that new vectors are now selected 
at intervals of P — 2tt/VL. This is the natural resonance 
timescale for the disk. We run two simulations with this 
altered forcing, both starting at the fiducial 20 orbital 
periods and running out until 150 periods, one having 
2/o and the other 8/o as their forcing amplitudes. We 
label these runs 2f$l and 8/r^, where the I indicates the 
long term nature of the forcing function. 

Our final two large s cale forcing simulati o ns are mo- 
tivated by the work of iLesur fc Longarettil (|2007| ) and 
iFromang et al.l (|2007f ) who demonstrate a link between 
the value of the magnetic Prandtl number v/rj and the 
saturation level of magnetic turbulence. We initialize 
two forced runs at 20 orbits with a 2/o forcing amplitude 
(as this seems to correspond to the demarcation between 
the standard MM behavior and that influenced by the 
hydrodynamic forcing) and run them to 150 orbits. For 
one run, labeled 2fop4 we increase the hyperviscous co- 
efficient by a factor of four to correspond to a four fold 
increase in the Prandtl number. For the other run, la- 
beled 2/op.25 we reduce the value of the hyperresistive 
coefficient to correspond to a four fold decrease in the 
Prandtl number. 

2.5. Small Scale Forcing 

In addition to examining the effects of large scale forc- 
ing we explore a (smaller) range of simulations where we 
input energy at a much smaller scale. Once again we set 
up a non forced MM simulations with conditions equiv- 
alent to those described above. We run this to 20 orbits 
and save the values at this point. We then allow the 
unforced MM to run to ~ 110 orbits and label this run 
Os/o where the s indicates the small scale forcing. 

We then repeat the process, described above, of de- 
termining the value of the forcing strength parameter 
s/o which produces a saturated kinetic energy level in a 
purely hydrodynamical run which is approximately equal 
to that produced in the fiducial MM run, except now 
forcing at k — 40. This value is chosen as it represents 
small scales within our inertial range which are not yet 
affected by dissipation due to numerical effects. Having 
obtained an approximate match between the kinetic en- 
ergy densities in the MRI-only and hydrodynamic runs, 
we label the forcing amplitude as s/o and implement a 
series of forced runs starting from the fiducial 20 orbit 
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Fig. 1. — Time and volume averages (equations l9l and l8)l of the 
transport coefficient, a (equation[TJ for the fiducial MRI run (0/o). 

point using multiples of this value. These runs are la- 
beled ls/o,2s/o, and 4s/o- 

3. RESULTS 

Due to the fluctuating nature of turbulent quan- 
tities both spatial and temporal averaging is needed 
in order to obtain reproducible values for quanti- 
ties such as the saturated ma gnetic field strength 
dWinters. Balbus fc Hawlevl l2003HPapaloizou fc Nelson! 
2003) . We present results in terms of the volume average 
as a function of time, given by 



< F(t) > = 



V 



F(t)dV, 



(8) 



< 



(9) 



and in terms of a running time average given by 

1 /■' 

F >st= — — / < F{t) > alt. 

We use the latter quantity to assess whether our runs 
are long enough to attain convergence. We note that 
our run length has been chosen to be just sufficient to 
yield clear discrimination between the different forcing 
cases — much longer runs would be required to measure 
quantities to percent level precision. 

3.1. Large Scale Forcing Runs 

3.1.1. Fiducial Run (0/o) 

We compare all of our runs to a fiduical non-forced 
MM configuration. This run displays the characteristic 
behavior of MM simul ations (|Brandenburg et "all 119951 : 
IBalbus fc Hawlevl [1991 . There is an initial exponen- 
tial growth phase which amplifies all components of the 
magnetic and velocity fields, which reach peak ampli- 
tude at ~ 3 orbits. After the peak there is a gradual 
reduction in all quantities until they saturate at a (fluc- 
tuating) level several orders of magnitude above their 
initial levels (except for B 2 Z which, in these zero-net flux 
simulations, saturates at a level below its initial value). 
The saturated state is sustained and all quantities re- 
main approximately constant until we stop the run at 
orbit ~ 150. 

Figure [1] shows the total transport coefficient a for the 
MM only case. The behavior is as described above: an 
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TABLE 1 

Model Parameters. Notes: ("'The maximal value of = P ga s/Pmagnetic, at t = 0; ( 6 'A C = 2nV az /V3f2; ( c )denotes a time and 

VOLUME AVERAGE FROM THE 20TH ORBIT UNTIL THE SIMULATION'S CONCLUSION 



exponential increase in a followed by a decay to reach 
a saturated state. The saturated value of a is around 
1.3 x 10~ 3 . As has been found in all prior simulations, the 
level of transport associated with the magnetic stresses is 
almost an order of magnitude higher than that associated 
with the fluid stresses. 

The MRI-only run serves two purposes. First, we use 
it as a reference for comparison with the forced runs. For 
this comparison, we compute running averages from the 
MRI-only run starting at orbit 20, thereby excluding the 
initial (unphysical) exponential growth and decay phase. 
Second, we use the snapshot generated at orbit 20 as an 
initial condition for the forced runs. These forced runs 
thereby model the effect of additional turbulence on an 
already saturated disk. We do not study how strong 
forcing would impact the linear growth phase of the MRI, 
as it seems unlikely that this is relevant to any physical 
system. 

3.1.2. Hydro Run (HL) 

We next ran several forced, non-MHD simulations, in 
order to find the value of the forcing amplitude (defined 
through /o) that yields hydrodynamic turbulence of the 
same strength as that obtained in the MRI-only run. 
MHD and non-MHD turbulence can be physically dif- 
ferent (jSridhar fc Goldre ich 1994J), so there is no unique 
definition of equivalent strength. We define an equiva- 
lent hydrodynamic run as one which yields close to the 
same value of the saturated kinetic energy density as the 
MRI-only run 3 . Numerically, we found that a value of 
fo = 0.0035 yielded reasonable agreement (at a level of 
about 10% between orbits 80 and 150). For this matched 
case Figure [2] shows how the running time averages of the 
kinetic energy compare in the MRI-only and hydro-only 
runs. Convergence of the kinetic energy in the box is 
rather slow in the hydrodynamic case (possibly due to 
the low level of viscosity, since lower resolution test runs 
converged more rapidly), but the level of the agreement 
at the end of the runs is sufficient for our purposes. 

Figure [3] shows the running time average of the fluid a 
in the hydrodynamic only reference run with this value 

3 Alternatively, we could define equivalence via the rate of energy 
input. This, however, is harder to measure in the MRI case. 
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Orbits (27T/0) 

Fig. 2. — Running time averages (equation |9j of the kinetic en- 
ergy for the unforced MRI run (0/o) and the hydrodynamic run 
(HL). In the hydrodynamic run forcing was used. The time aver- 
ages for the MRI run are begun at 20 orbits to prevent contam- 
ination from the initial exponential phase of the instability. The 
mean value of the time averages agree to within ~ 10% after orbit 
80. 

of /o. As expected for randomly forced hydrodynamic 
turbulence, the level of transport obtained is consistent 
with zero and certainly negligible compared with either 
the Maxwell or Reynolds stresses generated in the MRI- 
case. This result is important as we are explicitly looking 
at the effects of large scale hydrodynamic turbulence on 
the MRI. If a particular forcing function were to generate 
a non-zero value for a it would be difficult to separate 
its effects on transport from that of the MRI. 

3.1.3. Forced Runs (±/ , l/o, 2/ ,4/ Q , 8/ , 16/ ) 

Having determined the level for fo that matches the ki- 
netic energy between the MRI-only and hydro-only runs, 
we initiate six MHD runs with the amplitude of the forc- 
ing set at |/o, l/o, 2/o, 4/o, 8/0, and 16/ . We are 
particularly interested in the behavior of the magnetic 
fields and the angular momentum transport obtained in 
these simulations. 

Figure 0] shows the time evolution of the volume av- 
erage of B 2 Z for the fiducial unforced MRI simulation, 
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Fig. 3. — Time and volume averages (equations [9] and U of the 
transport coefficient, a (equation]!} for the purely hydrodynamical, 
large scale forcing simulation (HL). No transport is induced by the 
hydrodynamic forcing when it is employed on its own. 
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Fig. 4. — Evolution of B^, averaged over the simulation 
volume, for all MHD simulations including large scale forcing 
(0/o, l/o, 2/o, 4/ , 8/o, & 16/o). All forcing strengths above the l/ 
case induce noticeable increases in field strength. 

together with the results from the six simulations with 
variable forcing levels. The amplitude of the vertical 
magnetic field component is statistically similar (given 
the large fluctuations) for the 5/0 and l/o cases as com- 
pared to the unforced MRl-only run. For stronger forcing 
(2/o and above) B\ clearly exceeds that which would be 
generated by the action of the MR1 alone. Very roughly 
the energy in the vertical field scales linearly with the 
amplitude of forcing (recall that the forcing amplitude is 
defined via the energy in fluid motions). 

Our interpretation of the behavior of B z in the pres- 
ence of forcing is that the existence of external hydrody- 
namic turbulence introduces more turnover in the verti- 
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Fig. 5. — Evolution of B 2 /2/io, the total energy in the mag- 
netic field within the simulation volume. Results are shown for 
the unforced and all forced MHD runs in the large scale forcing 
simulations (0/o, l/oi 2/o, 4/o, 8/0, & 16/o)- For forcing amplitudes 
greater than 4/o we see a sharp rise in the strength of the magnetic 
field consistent with the idea that a new MRI exponential growth 
phase is initiated, subsequent to which the new saturation level is 
maintained. 

cal direction via random fluid motions. This increase in 
the "kneading" of the z component of the magnetic field 
by the forced turbulence increases B\. Since the genera- 
tion of vertical field by the MRI is relatively weak (this 
is reflected in the fact that, at saturation, only a small 
fraction of the magnetic energy in MRI simulations is 
in the z component) it is no surprise that modest levels 
of isotropic hydrodynamic turbulence suffice to generate 
higher B z . 

Figure [5] displays the total magnetic field energy den- 
sity for all the MHD simulations. For all runs we see an 
increase in the total field energy as /o increases. For 
the weaker forcing simulations (those up to 2/o) the 
timescale over which the disk adjusts to a new quasi- 
steady state (with a higher saturation level) is relatively 
long - typically tens of orbital periods. Somewhat differ- 
ent behavior is seen in the very strongly forced runs (4/o 
and above). In these cases the field adjusts almost im- 
mediately to a higher value once the forcing is turned on, 
and these stronger fields are then sustained on average 
for the duration of the simulation. What may be happen- 
ing here is that the enhanced B z created by the forced 
turbulence has "resparked" the rapid growth phase of 
the MRI allowing for more rapid adjustment of the field 
configuration. 

The presence of stronger magnetic fields in the stirred 
runs raises the question of whether the resulting turbu- 
lence (i) retains the character of the MRI, (ii) is more 
similar to non-magnetic turbulence, or (iii) is distinct 
from cither. We have attempted to quantify the charac- 
ter of the turbulence in two ways. First, we study the 
global characteristics of the turbulence, concentrating on 
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Fig. 6. — Running time average of a for all large scale forcing 
MHD simulations (0/o, l/o, 2/o, 4/ , 8/ , & 16/ ). 
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Fig. 7. — Running time average of the Maxwell and Reynolds 
components of a as well as the total a for the unforced MRI sim- 
ulation and the run with the strongest large scale forcing (16/o). 
The magnetic component of the transport remains dominant in 
both the unforced MRI run and in the most strongly forced case. 

the strength of angular momentum transport and how 
that stress is partitioned between Maxwell and Reynolds 
components. A notable feature of the MRI is that the 
stress is dominated by the Maxwell contribution, al- 
though the Reynolds stress is also non-zero. Second, we 
look at the structure of the turbulent flow as measured 
by the power-spectrum of various quantities. 



Fig. 8. — Power spectra of the kinetic energy (left panel) and the 
magnetic field energy density (right panel) in the large scale forcing 
simulations. The spectra have been averaged over orbits 130-140. 
HL is the hydrodynamic only large scale forcing simulation, 0/o is 
the unforced MRI simulation and the remaining curves show the 
forced, large scale simulations. The expected Kolmogorov slope 
has been plotted for the kinetic energy spectra and the location of 
the critical wavelength has been plotted in the panel corresponding 
to the magnetic field spectra. In all cases the total energy in the 
system increases with forcing amplitude. There is a transition be- 
tween MRI-like and hydrodynamic behavior visible in the kinetic 
energy spectra at large spatial scales. 

Figure [6] shows the time evolution of the equivalent a 
parameter for all of the runs that incorporate large-scale 
forcing. The division between magnetic and hydrody- 
namic stress is shown in Figure [7] for the unforced (0/o) 
and most strongly forced (16/ ) cases. Forcing boosts 
the strength of angular momentum transport once it is 
as strong (l/o) or stronger than the MRI in terms of the 
saturated equivalent kinetic energy. The derived a in- 
creases steadily with the forcing level up to the 16/o run, 
which displays a value of a about 6 times larger than 
the unforced MRI-only simulation. The fraction of the 
transport contributed by Maxwell stress as compared to 
Reynolds stress remains surprisingly constant between 
the unforced and maximally forced cases. Even though 
we are imposing external forcing - in some cases at quite 
a vigorous level - the transport is still determined by the 
magnetic field, and appears to retain the characteristic 
behavior expected from the MRI. 

We can also quantify the structure of the turbulence 
via power spectra of the velocity and magnetic fields. 
Figure [8] shows power spectra for both the kinetic en- 
ergy and (where appropriate) magnetic energy density 
for the hydrodynamic reference simulation, for the un- 
forced MRI run, and for all of the large-scale forcing 
simulations. For the hydrodynamic simulation (shown 
in the left panel Figure [5]) , we see that the kinetic en- 
ergy is dominated by the wavelength which corresponds 
roughly to the box size at which power is injected. For 
wavenumbers between k = 3 and k = 30 the spectrum is 
approximately of the Kolmogorov form (oc fc -5 / 3 ). For 
the MRI-only run the spectra of the kinetic energy and 
magnetic fields are quite similar - they are essentially 
flat at large scales (k < 20) with a subsequent sharp 
decline toward the dissipation scale. When forcing is in- 
cluded in the MRI simulation we see the same pattern 
of behavior discussed above in the context of a. For 
the relatively weakly forced runs (0.5/o and l/o) both 
the kinetic and magnetic power spectra are only slightly 
modified from the MRI-only case. For stronger forcing 
we find that the large-scale hydrodynamic structure of 
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the turbulence (quantified by the kinetic energy power 
spectrum) increasingly resembles that of non-magnetized 
turbulence. For the 16/o run the large-scale slope is rea- 
sonably close to the Kolmogorov expectation. The power 
spectrum of the magnetic fields, on the other hand, looks 
like a scaled version of the MRI-only run even up to the 
strongest forcing levels. There is some evidence for a 
shift of the peak power toward lower k, which would be 
consistent with the MRI operating at larger scales in the 
presence of stronger magnetic fields. We interpret these 
results as suggesting that the magnetic character of the 
turbulence remains similar (at all scales) to that gener- 
ated by the MRI even in the presence of strong hydro- 
dynamic forcing. The stronger magnetic fields generated 
by adding forcing then merely scale up the angular mo- 
mentum transport efficiency from the value expected in 
the absence of such effects. 

Figure [9] shows a slice of the density profile in the x-y 
(r-</>) plane at 150 orbits for simulations HL, 0/o, 4/o,and 
16/o- The dominant effect of the shear is clearly visible 
- the turbulence is stretched out into ribbons elongated 
in the azimuthal direction. Comparing the plots, about 
the only feature that is visually apparent is the increase 
in the large-scale power (here in the density field) as the 
strength of hydrodynamic forcing is increased. The most 
strongly forced run displays substantially greater den- 
sity contrasts, and it is clear that the large-scale mor- 
phology of the flow has been altered significantly from 
the MRI-only case. Nevertheless, as discussed above, on 
the (somewhat smaller) scales that are most important 
for angular momentum transport the MHD properties of 
the flow are apparently not significantly modified. 

To quantify the relationship between the strength of 
imposed forcing and the resulting physical quantities (a, 
kinetic energy density etc) we measure the ratio of the 
mean quantities in the forced to the unforced run. This 
ratio, defined as, 



T(fo) = 



< F(t&nah fo) >8t 



(10) 



< F(t Rnah f = 0) > st 

is plotted in Figure \W\ for a, its decomposition into 
Maxwell and Reynolds stresses, the total magnetic field 
energy, and the total kinetic energy. It is evident that, 
to a good approximation, the additional energy injected 
via the forcing simply raises the saturated value of the 
kinetic energy in the same manner as one would expect 
for purely hydrodynamic runs - i.e. the I6/0 run, which 
would have 16 times the kinetic energy of the MRI-only 
run in the hydrodynamic limit, does indeed yield that en- 
hancement in the presence of MHD. The fact that these 
are MHD runs, in which the dissipation properties of 
turbulence could in principle differ significantly from the 
hydro only case, does not appear to make a substan- 
tial difference. This is similar to the situation studied 
in the context of molecular clouds — where MHD and 
hydrodynamic turbulence have comparable decay rates 
(|Mac Low et alj|1998t fStone. Ostriker fc Gammielll998h 
— though here the fluid motions are highly subsonic. 

As noted previously, all of the magnetic quantities, to- 
gether with the Reynolds stress (which is an approxi- 
mately fixed fraction of the total stress) increase in uni- 
son as the forcing strength increases. The relationship 
is roughly linear in the value of the forcing parameter 
fo, though with a smaller slope than the 1:1 relation- 



ship seen for the kinetic energy. A forcing strength that 
yields an order of magnitude increase in the saturated ki- 
netic energy density results in about a factor 4 increase 
in the the magnetic field energy density and in all of the 
components of a. 

3.1.4. Large Scale Runs With Different Prandtl numbers 

A number of recent studies provide evidence that the 
magnetic Prandtl number, 



Pm = 



V 



(11) 



is an important parameter whose value affects the sat- 
uration level of MHD turbulence initiated by t he MRI 
(|Lesur &; Longarettil [2007t iFromang et al.ll2007t ) 4 . This 
dependence is of interest since the predicted value of Pm 
varies substantially between the inner regions of black 
hole accretion disks (where Pm 3> 1) and the outer re- 
gion (|Balbus fc Henril r2008). In all of the simulations 
performed above our hyperviscous Prandtl number i>h / rjh 
was equal to one. To explore the affects of varying the 
the Prandtl number we ran two simulations, one in which 
we increased the value of our viscosity coefficient by a 
factor of four and one in which we reduced our resistive 
coefficient by a factor of four. This change in coefficients 
allowed for an effective span of 16 in Prandtl number. We 
ran both of these simulations with a forcing amplitude of 
2/o as this marked the lowest forcing level for which sig- 
nificant differences with the unforced runs were observed. 
We ran these simulations from our fiducial start point at 
20 orbits out to 150 orbits and labeled them 2/ p4 and 
2/op.25. 

Figure 1111 shows the volume averaged magnetic en- 
ergy density for simulations 2/ , 2/ p4, and 2/ p.25. 
In this comparison our 'fiducial' model is the 2/o run 
with a Prandtl number of one. We observe sustained 
turbulence out to 150 orbits in all cases, implying that 
our effective Reynolds number (which should be boosted 
somewhat by our use of hyperviscosity rather than a 
regular viscosity) is high enough that sustained turbu- 
lence is possible even for Prnj, = 0.25, i.e. that we 
are in the larg e Reyn olds number regime described by 
IFromang et al.l (|2007f K In fact, no clearly significant dif- 
ferences are seen between the Prnj, = 0.25 and Pmi, = 1 
runs. For Pnih = 4, on the other hand, we see an increase 
of about a factor of 2-2.5 in the magnetic field energy and 
angular momentum transport efficiency. Although the 
use of hyperviscosity and hyperresistivity means that our 
runs are not equivalent (and i ndeed are l ess physical) to 
those of iLesur fc LongarettH (|2007| ) and IFromang et al.1 
(|2007l ) this trend is consistent with those previous stud- 
ies which show an increase in saturation level with in- 
creased Prandtl number. We borrow the same interpre- 
tation: values of Pm > 1 result in strong viscous damp- 
ing at scales larger than the magnetic field dissipation 
scale. This suppresses magnetic field dissipation at small 
scales, and favors an inverse cascade of magnetic energy 
to larger scales and an attendant increase in a. 

While we see an increase in the saturation level for 
the higher Prandtl number simulation the the underly- 
ing behavior remains the same as that seen in the large 

4 The value of Pm als o affects the properties o f non-rotating, 
drive n MHD turbulence ( Schckochihin ct al. 2004; Iskako v~et al.l 
120071) . 
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Fig. 9. — Density slices at 150 orbits in the xy plane taken from simulations HL (upper left), 0/o, 4/rj, & 16/rj. 
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Fig. 10. — The average values of the magnetic field energy den- 
sity, kinetic energy and components of a from the large scale forc- 
ing runs, normalized to the results of the unforced MRI simulation. 
The forcing strength is represented on the x-axis by the value of 
the normalization parameter for the forcing /o. 

scale forcing runs. For both runs with altered Prandtl 
numbers the transport is still dominated by the Maxwell 



Fig. 11. — Evolution of the magnetic field energy density 
(B 2 /2/^o) for the simulations with altered Prandtl numbers (u/r]). 
All three runs used 2/o as the forcing amplitude. 

stresses and the proportion of Maxwell to Reynolds stress 
remains the same. We find no significant differences in 
the behavior of the field components or the power spectra 
except that they are amplified in the the higher Prandtl 
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Fig. 12. — Evolution of the magnetic field energy density 
(B 2 /2/io) f° r the simulations with orbital period forcing variations. 
The results are compared to those from the equivalent large scale 
simulations in which the forcing was varied on the timescale of the 
simulation timestep. The results marked with an '1' denote the 
orbital period simulations. No significant differences are seen. 

number run. 

3.1.5. Large Scale Runs With Orbital Period Forcing 

In the large-scale forcing runs described above the vec- 
tors which describe the instantaneous turbulent driv- 
ing are chosen anew at each timestep. The forcing is 
then temporally uncorrelated with any characteristic fre- 
quency of the disk. To check whether the results depend 
upon this feature we repeated some of the runs with the 
correlation time for the forcing set to be the orbital pe- 
riod 5 . We ran simulations at two times and eight times 
the reference forcing amplitude with orbital period forc- 
ing (2/nJ and 8/ ?) and compared them to the standard 
runs with time step forcing variations. Figure fT2l shows 
the volume averaged magnetic energy density in the sim- 
ulations with orbital period forcing variation compared 
to their standard counterparts. We find no discernible 
differences in the magnetic field, transport coefficients, 
or spectra in the case of orbital period variations. Both 
the behavior and the saturation amplitudes agree. Once 
again the Maxwell stresses dominate the transport. 

3.2. Small Scale Forcing Runs 

We have also considered the effects of forcing on scales 
A <C H, Once again we begin by finding the strength 
of forcing that yields a saturated kinetic energy density 
that matches that obtained in a fiducial MRI-only simu- 
lation, except now we use a forcing function that injects 
power at wavenumbers around k — 40. The match we 
obtained was good at the » 20% level. We found that 
the residual hydrodynamic a in the purely hydro run 
was substantially higher than in the case with large scale 
forcing, but still negligible for our purposes (being two 
order of magnitude below the MRI value). 

Figure [T3] and figure [TJ] show the time evolution of B% 
and B 2 /2fiQ for the small scale forcing simulations. We 
considered a more limited range of forcing strengths be- 
tween ls/o and 4s/q. Within this range, there are some 
differences between the behavior of the disk subject to 
large and small scale forcing. Most notably, although 

5 We adopted this timescale since it is a well-defined character- 
istic frequency for an accretion disk. Since it is also the resonant 
frequency, any effects due to changing the timescale are likely to 
be maximized by adopting orbital period scale forcing. 



Fig. 13. — Evolution of B 2 for the simulations with small scale 
forcing (0s/ , ls/o, 2s/ , &4s/ ). 
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Fig. 14. — Evolution of the magnetic field energy den- 
sity B 2 /2fiQ for the simulations with small scale forcing 
(Os/ ,ls/o ! 2s/o ! &4s/o). 
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Fig. 15. — Running time average of a for the simulations with 
small scale forcing (Os/o, ls/o, 2s/o, &4s/o). The lowest two forc- 
ing strengths result in a suppression of angular momentum trans- 
port. 

relatively weak levels of forcing (the 2s/o run) do in- 
crease the value of the vertical magnetic field, this is 
not accompanied by any significant increase in the to- 
tal magnetic field. The strongest level (4s /o) of forcing 
does amplify all components of the magnetic field. We 
also observe larger fluctuations in magnetic field strength 
as compared to the equivalent large-scale forcing simula- 
tions. 

The stress within the simulated patch of disk largely 
follows the total magnetic field strength. FigurefTolshows 
the running time average of the equivalent a parameter 



The MRI in stirred disks 



11 




10 100,1 10 



Fig. 16. — Power spectra of the kinetic energy (left panel) 
and the magnetic field energy density (right panel) for simulations 
with small scale forcing. In these runs the forcing is applied at a 
wavenumber k = 40. The spectra have been averaged over orbits 
95-105. HS is the hydrodynamic only simulation with small scale 
forcing, Os/o is the unforced MRI run, while the remaining curves 
correspond to the small scale, forced MRI runs. Unlike in the large 
scale forcing runs Kolmogorov spectra are not observed at large 
scales, and there is a marked accumulation of power close to the 
injection wavenumber. Power in the large scale magnetic field is 
essentially unaltered by the forcing except in the highest amplitude 
forcing run. 

for the different simulations. The average a is lower in 
the weakly forced runs (ls/o and 2s/o) than in the un- 
forced simulation, and only shows significant enhance- 
ment in the 4s/o case. In this case it is clear that the 
value of a has not converged by the end of the simula- 
tion. We note that in all cases - including those where 
the forcing appeared to suppress the MRI - the transport 
remains dominated by Maxwell stresses. 

Figure [16] shows the power spectra of the kinetic energy 
(left panel) and magnetic energy density (right panel) for 
all of the small scale forcing simulations. The accumu- 
lation of power in the kinetic energy spectrum at scales 
close to the forcing scale of k = 40 is very obvious in all 
the forced runs, and on this scale the power in fluid mo- 
tions greatly exceeds that which would be generated by 
the MRI. In the hydro-only run power leaks out to larger 
scales, and results in a roughly flat power spectrum out 
to the box scale at k = 1. The spectra of the magnetic 
field energy are rather different. In the two most weakly 
forced runs the forcing increases the power close to the 
forcing scale, but diminishes or leaves roughly unaltered 
the power on larger scales. In particular, for these runs 
there is no increase in the power on the scales - com- 
parable to the most unstable linear MRI modes - that 
are presumably most important for field evolution in a 
disk. Only for the strongest level of forcing (4s/o) do we 
see a clear increase in the power across all of the scales 
accessible to the simulation. 

Based on these results, we conjecture that the impact 
of hydrodynamic forcing on the MRI varies depending 
on the effect that the forcing has on magnetic fields with 
scales comparable to the most unstable MRI wavelengths 
in the disk. Large scale forcing efficiently enhances the 
vertical magnetic field in the disk at quite moderate am- 
plitudes, and essentially always boosts the strength of 
angular momentum transport. Power from small scale 
forcing primarily affects the structure of the turbulence 
at scales close to and below the injection scale, and, at 
least if the forcing is weak, relatively little power cas- 
cades back to large scales. There is no enhancement 



of the field strength on the scales most important for 
growth of the MRI. Rather, the imposed turbulent mo- 
tion on small scales acts to suppress the MRI, perhaps 
by enhancing magnetic field dissipation (the inverse of 
the Prandtl number effect discussed above), or perhaps 
by directly destroying the correlation between the radial 
and azimuthal field that is responsible for the Maxwell 
stress. Only if the small scale forcing is rather strong 
does enough power reach larger scales, at which point 
the MRI can be enhanced via the same mechanisms as 
apply in the large scale forcing runs. 

4. DISCUSSION 

In this paper we have used local shearing-box simula- 
tions of accretion disks to study the interaction between 
the magnetorotational instability and hydrodynamic tur- 
bulence in disks where both coexist. We have studied 
both large scale forcing, in which energy is injected at 
scales A ~ H, the disk scale height, and to a more lim- 
ited extent small scale forcing where A <C H. We find 
that the effect of the additional energy input from the 
hydrodynamic turbulence on the MRI varies depending 
upon both the strength of the forcing and, to some ex- 
tent, on its spatial structure. 

For large scale forcing the results are straightforward. 
When hydrodynamic forcing is "weak" - in the sense that 
the hydrodynamic forcing, on its own, yields a saturated 
kinetic energy density in the disk that is less than or 
equal to that generated by the MRI alone - the MRI is 
essentially unaffected at the level of precision accessible 
to our simulations. This result is no surprise. The MRI is 
a robust instability, which is present within differentially 
rotating flo ws containing almost arbi trary magnetic field 
geometries (jBalbus fc Hawlevl Il998( ). Low amplitude 
random perturbations do not affect it. When stronger 
forcing is imposed, we find that both the saturation value 
of the magnetic field and the strength of angular momen- 
tum transport can be substantially boosted. The flow 
retains many of the characteristics of the MRI — such 
as a similar ratio between magnetic field strength and 
a, and a similar ratio of Maxwell to Reynolds stresses, 
even in a regime when the hydrodynamic forcing (which, 
by itself, yields no transport at all) is formally domi- 
nant. The application of heuristic dynamo ar guments 
to the MRI is suspect (jBalbus fc Hawleylll99c? ). but we 
tentatively attribute the numerical behavior to the more 
efficient regeneration of vertical field in the presence of 
hydrodynamic turbulence. 

The small scale forcing results are more nuanced, and 
given the limited number of simulations we have per- 
formed should be regarded as preliminary. In this regime 
we find that only the strongest level of forcing boosts the 
strength of angular momentum transport, whereas lower 
forcing levels actually suppress transport. We attribute 
this different behavior to the fact that there are two im- 
portant considerations that affect the saturation ampli- 
tude of the MRI. One is the strength of the vertical mag- 
netic field on relatively large scales (typically a fraction 
of H), similar to the most unstable linear MRI wave- 
lengths. Enhancement of the vertical field on this scale 
- which is readily accomplished with large scale forcing 
but which requires an inverse cascade in the small scale 
case - boosts the strength of the MRI. The second is 
the dissipation scale. Turbulent driving that increases 
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the amplitude of kinetic energy at or near this scale may 
enhance magnetic field dissipation, ultimately reducing 
the saturation level of magnetic fields in the disk as a 
whole. This is the inverse of the physical process invoked 
to explai n the dependence of a on the magnetic Prandt l 
number (jLesur & Longaretti 2007; F romang et alJl2007[ ). 
Small scale forcing may also directly destroy the correla- 
tions between B r and that result in Maxwell stress. 

Our results allow us to infer which additional physical 
effects are likely to be able to affect the MRI in accretion 
disks. We first observe that if the additional turbulence 
is "powered" by the MRI itself (i.e. if the turbulence 
derives energy from the gravitational potential well as 
a result of MRI-driven angular momentum transport), 
then generically it is unlikely to be as powerful as the 
MRI. Such turbulence will have at most a small effect 
on the magnitude and character of angular momentum 
transport. As an explicit example, we would not ex- 
pect an MRI-active disk that was additionally unstable 
to convection to differ much from one in which the ver- 
tical structure was stable against convective motions. 

We can also consider sources of turbulence that are in- 
dependent of the MRI, in the sense that they would exist 
even in a (hypothetical) disk that was absent magnetic 
fields entirely. Physical examples include self-gravity, 
thermal instability, and turbulence stirred up by the 
interaction between the gaseous and solid components 
of protoplanetary disks. There is no reason why these 
sources of turbulence should not generate fluid motions 
of greater amplitude than those produced by the MRI 



(this is likely to be the case in self-gravitating disks near 
the fragmentation threshold, and locally in some regions 
of protoplanetary disks). In this regime, our results sug- 
gest that the MRI will have a non-trivial interaction with 
the hydrodynamic turbulence. If the forcing occurs at 
scales comparable to H or larger, we find that the inter- 
action will likely boost the strength of angular momen- 
tum transport. Even for quite violent forcing - up to an 
order of magnitude in excess of that required to produce 
parity with the MRI - we find that it is most accurate to 
think of the coupled system as one with boosted MHD 
turbulence, rather than as a hydrodynamic system pas- 
sively advecting magnetic fields. Conversely, small scale 
forcing, unless it is very strong, may actually suppress 
the saturation level of magnetic fields in the disk and 
their associated angular momentum transport. 
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